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TECHNICAL PAPER 


FURTHER DEVELOPMENTS IN EXACT STATE RECONSTRUCTION IN 
DETERMINISTIC DIGITAL CONTROL SYSTEMS 

I. INTRODUCTION 


Books on modem digital control systems usually address the problem of controlling a continuous- 
time plant driven by a zero-order-hold with a sampled output as shown in Figure 1 (for example, see 
Reference 1 , p. 126). A common solution to this problem is to reconstruct the state of the system at the 
sampling instant using a state observer and then feed back the reconstructed state [1, p. 195]. However, 
the state observer has two undesirable characteristics. First of all, it is a dynamical system in itself and, 
hence, adds additional states and eigenvalues to the system, which can affect system stability. Second, as 
a consequence, the reconstructed state is normally an approximation to the true state and is usually not a 
good one early in the state reconstruction process. Recently, Polites developed a new approach to state 
reconstruction which has neither of these problems [2]. Subsequently, he extended this work and 
developed what he called the ideal state reconstructor [3], It was so named because: if the plant 
parameters are known exactly, its output will exactly equal the true state of the plant, not just approxi- 
mate it. Besides that, it adds no new states or eigenvalues to the system. Nor does it affect the plant 
equation for the system in any way; it affects the measurement equation only. It is characterized by the 
fact that measurements prefiltered by a multi-input/multi-output moving-average (MA) process [4] are 
utilized in the state reconstruction process. Now, in this paper, a more-general version of the ideal state 
reconstructor is presented. It allows standard instantaneous measurements to supplement the MA- 
prefiltered ones in the state reconstruction process. Useful in the development of this more-general ideal 
state reconstructor are some results to date for continuous-time plants driven by a zero-order-hold with 
sampled outputs. These are reviewed in Section II, prior to the development of the more-general ideal 
state reconstructor presented in Section III. Section IV presents an example which illustrates the proce- 
dure for choosing the parameters in it. Section V contains the conclusions and recommendations. 
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Figure 1 . Continuous-time plant driven by a zero-order-hold 
with instantaneous measurements. 




II. PRELIMINARY 


For the plant in Figure 1 , x(t) is an nx 1 state vector, u(kT) is an rx 1 control input vector, yi(kT) is 
an mxl output or measurement vector, F is an nxn system matrix, G is an nxr control matrix, and Q is an 
mxn output matrix. Since y^k) = C T x(k), where k is the usual shorthand notation for time kT, y ( (k) 
represents an instantaneous measure of the system at the sampling instant kT. Hence, the plant in Figure 
1 can be regarded as having instantaneous measurements for outputs. It is well known that this system 
can be modeled at the sampling instants by the discrete state equations [1, p. 126] 


x(k+ 1) = 

Ax(k) + Bu(k) 

(1) 

yi(k) = c, 

x(k) , 
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where 


<t>(0 [(sI-F)- 1 ] , 

(3) 

A = 4>(T) , 
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and 


B = 


r T 

o 


G 


(5) 


4>(t) is the nxn state transition matrix. A and B are the nxn system matrix and the nxr control matrix, 
respectively, for the discrete state equations (1) and (2). 

A and B can be determined analytically using equations (3) through (5). An alternative approach, 

t 

which is also quite suitable for numerical computation, is as follows [5]: <\>(t) and f <j>(\)dA. can be 

o 

expressed in the form of matrix exponential series as 
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and 
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<t>(\)d\ = 2 
i = 0 


F t M 

7T+T)! 
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respectively. From equations (6) and (7), 


<|>(t) = I + F 


J^/W)dX 


( 8 ) 


where I is an nxn identity matrix. Hence, f 4>(X)d\ can be determined using equation (7) with t = T 

o 

and this result substituted into equation (8) to get 4>(T). With these results, A and B can be found using 
equations (4) and (5). 

Now consider the plant in Figure 2, which is a generalization of the one in Figure 1 . In addition to 
the instantaneous measurement vector yi(kT), the plant in Figure 2 has the measurement y F '(kT) gener- 
ated as follows. First, the continuous-time output z(t) is sampled every T/N seconds. Every N samples 
are multiplied by the weighting matrices Hj, j = 0,1,..., N-l, and then summed to generate the output 
y F (kT), every T seconds. Functionally, this is equivalent to passing the discrete measurements generated 
every T/N seconds through a multi-input/multi-output MA process with coefficient matrices Hj, j = 
0,1,..., N-l [4]. The output of the MA prefilter is sampled every T seconds to generate y F (kT). Then, 
y F (kT) has subtracted from it E_u[(k- 1 )T] , where E. is a constant matrix, to produce the modified MA- 
prefiltered measurement vector y F '(kT). Finally, y F '(kT) is catenated with y,(kT) to form the total mea- 
surement vector yx(kT). In Figure 2, C F is a pxn output matrix and z(t) is a pxl vector. The weighting 
matrices Hj, j = 0,1,... ,N-1 are each qxp. Hence, y F (kT) and y F '(kT) are qxl vectors. Since yi(kT) is an 
mxl vector, it follows that y-r(kT) is an (m + q)xl vector. Since u[(k-l)T] is an rxl delayed input vector, 
E. is a qxr matrix. 

Previously, Polites [6] showed that when 


E. = HP , (9) 

where H is a qx(Np) matrix given by 

H = [Ho J H, | ... | H n .,] (10) 
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and 3 is the (Np)xr matrix 



C F 

4>(X)dx] G 

P = 

-(T/N) 

C F [J 4>(\)dx] G 
0 


• 

• 

C F | 

r -(N-l)(T/N) 

J 4>(X)dxJ G 
0 J 


the discrete state equations for the plant in Figure 2 become 

x(k+l) = Ax(k) + Bu(k) (12) 


y-r(k) = 


’yi(k) " 


"c,“ 

yp'(k) 

m m 


D_ 


?(k)=c T x(k) , 


(13) 


where D. is a qxn matrix given by 


D. = Ha 


(14) 


and a is the (Np)xn matrix 


a 


" C F <j>(0) 

Cp ?<- N-> 

. Cf * [-(N-l) i] ■ 


(15) 
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From equation (13), 


Cy 


C, 


D. 


( 16 ) 


where C T is an (m + q)xn matrix. 

E. and D. can be evaluated analytically using equations (3), (9) through (11), (14), and (15). An 
alternative approach, which can be either analytical or numerical, is as follows. Let t = -j(T/N), where j 

-j(T/N) 

= 0,1 ,...,N-1, and use equation (7) to determine J <J>(\)d\,j = 0,1 ,...,N-1 . Substitute these results 

o 

into equation (8) to get <f>[-j(T/N)], j = 0,1,...,N-1. At this point, E. and D.can be found using equations 
(9) through (11), (14), and (15). 

III. THE MORE-GENERAL IDEAL STATE RECONSTRUCTOR 


A general block diagram of the plant and the more-general ideal state reconstructor is shown in 
Figure 3. Notice the similarity between Figures 2 and 3. By virtue of this, if E. is given by equation (9), 
then equations (12) to (15) define the discrete state equations for the system in Figure 3 up to the output 
y T (k). Proceeding further, y T '(k) is related to y T (k) by the expression 


y-r'00 — (Ct T Cj)" 1 Cj y-r(k) 


(17) 


However, for equation (17) to be meaningful, (Cx T C T )'' must exist, and this occurs only when (C T T C T ) 
is nonsingular. Recall that C T is an (m-l-q)xn matrix. If (m + q) 3* n and C T has maximal rank (i.e. rank 
n), then (Cx T C T ) is positive definite and therefore nonsingular [7]. Hence, equation (17) requires that 
(m + q) 3= n and rank (C T ) = n for it and the more-general ideal state reconstructor to be meaningful. 
Assuming this is the case, it follows from equations (12), (13), and (17) that the discrete state equations 
for the system in Figure 3 are 


x(k+l) = Ax(k) + Bu(k) (18) 

y T '(k) = x(k) . (19) 


Hence, the output of the more-general ideal state reconstructor, y T '(kT), exactly equals the true state of 
the plant, x(kT). Consequently, if one is given the plant in Figure 1 and modifies it to conform to Figure 
3, he can exactly reconstruct the state of the plant without adding any new states, eigenvalues, or 
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dynamics to it, since the plant equation (18) for the system in Figure 3 is identical to the plant equation 
( 1 ) for the plant in Figure 1 . In Figure 3 , exact state reconstruction is achieved when E_ is given by equa- 
tion (9) and Ct is given by equation (16) where D. is given by equation (14). Of course, C T must satisfy 
the requirements previously imposed on it. 

One of these is that the dimensions of Cj, namely (m + q)xn, satisfy the relationship (m + q) 3= n. 
This can also be written as q 3= (n-m). Hence, the number of rows, q, in the weighting matrices Hj, j = 
0,1 ,...,N-1 must equal or exceed the number of states, n, in the plant state vector, x(t) or x(kT), minus 
the number of elements, m, in the instantaneous measurement vector yi(kT). Since q can be chosen 
arbitrarily, this requirement can be readily satisfied. The other requirement on C T is that it have rank n. 
One approach to satisfying this requirement is as follows. Recall that C T is defined by equation (16) 
where D. is defined by equation (14). Assuming Cj is given, then one can choose D. so that C T has rank n 
and then find H to give the desired D_. One solution to the problem of finding H to give the desired D_ is 
to let 


H = D. (a T a )' 1 a T . 


( 20 ) 


This follows from equation (14). However, like before, this requires that (a T a) be nonsingular. Recall 
that a is an (Np)xn matrix. If (Np) 3 n, or equivalently N 3= n/p, and a has maximal rank (i.e., rank n), 
then (a T a) is nonsingular. The first requirement can be easily satisfied because the number of weighting 
matrices N, where the weighting matrices are Hj, j = 0,1 ,. . . ,N-1 , can be arbitrarily chosen so that N 3 
n/p. Recall that n is the number of states in the plant state vector, x(t), and p is the number of elements in 
the output vector z(t). 

In summary, the procedure to achieve exact state reconstruction with the more-general ideal state 
reconstructor is as follows. Given the plant in Figure 1, modify it to conform to Figure 3. Choose the 
number of rows, q, in the weighting matrices Hj, j = 0,1,...,N-1 so that q 3= (n-m) where n is the 
number of states in the plant state vector, x(t) or x(kT), and m is the number of elements in the instantan- 
eous measurement vector y r (kT). Choose the number of weighting matrices, N, so that N 3* n/p where p 
is the number of elements in the output vector z(t). Assuming the (Np)xn matrix a, defined by equation 
(15), has maximal rank (i.e., rank n), let H be given by equation (20) where D. is chosen so that the 
(m-l-q)xn matrix Ct, defined by equation (16), has maximal rank (i.e., rank n). The weighting matrices 
Hj, j = 0,1,...,N-1 are found by partitioning H as in equation (10). Finally, let E. in Figure 3 be given by 
equation (9) where 3 is defined by equation (11). The discrete state equations for the system in Figure 3 
are now given by equations (18) and (19). Hence, the output of the more-general ideal state reconstruc- 
tor, y T '(kT), exactly equals the true state of the system, x(kT). 

In the event there are no instantaneous measurements, then y t (kT) is a null vector, Q is a null 
matrix, and the more-general ideal state reconstructor in Figure 3 degenerates to the one presented in 
Reference 3. In this case, the methods described in Reference 3 for choosing the parameters in the ideal 
state reconstructor to achieve exact state reconstruction, as well as the one presented here, are applicable. 
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IV. AN EXAMPLE 


Consider the double integrator plant driven by a zero-order-hold as shown in Figure 4. Manipulat- 
ing the plant in Figure 4 into the state variable format of Figure 1 yields 


F 



1 

0 


G 



( 21 ) 


( 22 ) 


and 


C, =[10]. 


(23) 


U( kT ) ^ 

ZOH 

u(t) = x 2 


*2=*1 


X, 

/ y,(kT) 




T 


Figure 4. Plant for the example. 


Since F, G, and C T are nxn, nxr, and mxn matrices, respectively, it follows from equations (21) to (23) 
that n = 2, r = 1, and m = 1. Using equations (21) and (22) and the formulas presented in Section II, 


4>(t) = 


1 • t 
0 1 
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r 




t 


> 


(24) 


(25) 
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A = 


1 

0 


T 

1 



To exactly reconstruct the state of the plant in Figure 1 with the more-general ideal state recon- 
structor, modify it to conform to Figure 5. Comparing Figures 3 and 5, it is apparent that 


C F = [1 0] 


(26) 


Since C F is a pxn output matrix, it follows that p = 1 . The requirement q 5* (n-m) can be satisfied by 
letting q = 1 • The requirement N 2* n/p can be satisfied by letting N = 4. Now a and (3 can be evaluated 
using equations (11), (15), (22), (24), (25), and (26), and are found to be 
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a = 



and 
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T 2 



4T 2 

32 

9T 2 

32~ 


(27) 


( 28 ) 
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respectively. In equation (27), eliminating any two rows forms a 2x2 matrix with nonzero determinant, 
assuming of course that T > 0. Hence, rank (a) = 2 = n and so (a T a) is nonsingular. Consequently, 
(a T a)"' a T exists and is found to be 


( |o ) (“>) (">) (■ '.) 
.(*) (w) ( 4 ) ( 4 ). 


(29) 


using equation (27). Since (a T a)"' a T exists, H can be given by equation (20) where D. needs to be 
chosen so C T has maximal rank. Recall that C T is an (m + q)xn matrix. In this example, m = q = 1 , n = 
2, and so C T is a 2x2 matrix. Since C T is defined by equation (16) where Q is the 1x2 matrix in equation 
(23) for this example, D_ must be a 1x2 matrix. If D_ is chosen to be 


D. = [0 1] 


(30) 


then 


1 


Cj 



0 

1 


and rank (Ct) = 2 = n since det(C T ) + 0. Hence, (C T T C T ) is positive definite and therefore nonsingular. 
Consequently, (Ct t Ct)"' C t t exists and is found to be 


(C-i^Ct)" 1 C t t 


1 0 
0 1 


From equations (10), (20), (29), and (30), 


H = [H 0 j H, | H 2 i H 3 ] = 


(4-) ! (4-) i ( if ) (■ if) 


(31) 
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which reveals the weighting matrices Hj, j = 0,1, 2, 3. From equations (9), (28), and (31), 


E. = - 3T/8 . 


The more-general ideal state reconstructor is now completely defined for this example. 


V. CONCLUSIONS AND RECOMMENDATIONS 


This paper has presented a more-general version of the ideal state reconstructor for deterministic 
digital control systems previously developed by Polites [3]. It is called an ideal state reconstructor 
because, unlike the popular state observer, it adds no new states, eigenvalues, or dynamics to the 
system, and, consequently, will not alter the stability of the system. In fact, adding the ideal state recon- 
structor to the system will not affect its plant equation in any way. It affects the measurement equation 
only. Also, if the plant parameters are known exactly, the reconstructed state will exactly equal the true 
state of the system, not just approximate it. The ideal state reconstructor is characterized by the fact that 
measurements prefiltered by a multi-input/multi-output moving-average (MA) process [4] are utilized in 
the state reconstruction process. It is called a more-general version because, unlike the original version 
[3], it allows for instantaneous measurements to supplement the MA-prefiltered ones in the state recon- 
struction process. The disadvantages of either version are that measurements must be made and calcu- 
lations must be performed more frequently than with the state observer. Fortunately, this is not the 
problem it was 20 years ago, considering the speed of today’s digital computers. 

If the research in this paper is extended, two approaches are recommended. One is to explore 
advanced methods for choosing the parameters in the ideal state reconstructor. For example, the weight- 
ing matrices might be selected so the MA prefilter acts as a multi-input/multi-output low-pass filter for 
the case where measurement noise is present. The other approach is to investigate the robustness of the 
ideal state reconstructor and see how it compares with the state observer’s. Specifically, the following 
questions should be addressed. What effect do modeling errors in the plant have on the ideal state recon- 
structor, and how does this compare with the state observer? What effect do plant process and mea- 
surement noise have on the ideal state reconstructor and how does this compare with the state observer, 
or even the Kalman filter? How can the robustness of the ideal state reconstructor be improved? Increas- 
ing the number of weighting matrices, N, may be one possibility. Catenating the ideal state reconstructor 
with a state observer, or a Kalman filter, may be another. This might produce a composite estimator 
which is better than either the ideal state reconstructor, the state observer, or the Kalman filter alone. 
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